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1. Introduction 

One of the most precise determinations of the strong coupling constant a s (Mz) relies on the next- to- 
next-to-leading order (NNLO) QCD analysis of the world deep -inelastic data [1,2]. Here currently 
the heavy flavor corrections to the structure function F2(x,Q 2 ) and transversity are known for the 
first Mellin moments N = 2, .., 10(14) [3,4]. ! In the asymptotic region Q 2 3> m 2 the heavy flavor 
Wilson coefficients can be represented in form of convolutions [7] of massive operator matrix 
elements (OMEs) and the massless Wilson coefficients [8]. In case of the charm quark contribution 
the corresponding region is given by Q 2 /m 2 > 10. To carry out complete NNLO QCD analyses in 
this region the heavy flavor Wilson coefficients have to be known for general values of N. This also 
applies to precision measurements of the charm quark mass using the world deep-inelastic data [9]. 

Since 2010 the systematic calculation of the asymptotic massive Wilson coefficients at gen- 
eral values of N have been carried out. Five massive Wilson coefficients contribute to F%{x, Q 2 ) 
at NNLO [3]. Furthermore, there are other massive OMEs needed to compute the matching coef- 
ficients in the variable flavor number scheme in which the heavy quarks are assumed to decouple 
singly [3, 10, 11]. There are yet other contributions at NNLO from graphs containing both a mas- 
sive charm and a bottom quark line, which extend the former representation in Ref. [3] and are 
necessary because of the fact that charm quarks do not yet become massless at the mass scale of 
bottom quarks since m 2 /m 2 b pa 1/10 only. Results on first moments for these contributions are ob- 
tained in [12-14]. The extension of the renormalization of the massive OMEs is given in [14]. 
The logarithmic contributions at general values of N for the contributions to F2(x,Q 2 ) are avail- 
able [15- 17]. 2 Furthermore, the NNLO heavy flavor Wilson coefficients in the asymptotic region 
were calculated for the structure function Fl(x, Q 2 ) [3,20]. 

Two of the five massive Wilson coefficients contributing to F%{x, Q 2 ) at NNLO are known 
completely [21]. They are of OfofTpCp^). Likewise, these contributions to the further three 
Wilson coefficients and transversity were calculated for these color coefficients in [21]. Also the 
complete contribution 0{T%C Ai f) for the OMEs A^,A^ s e ,A™ g are available [12]. 

In these proceedings we report on the calculation of the gluonic OMEs 0(nfTpCf^) [22] in 
Section 2 and on first results in 0(TpCfa) for this channel in Section 3. Furthermore, the scalar 3- 
loop integrals for all ladder type integrals were calculated [23,24]. An extension of the method [25] 
to calculate finite Feynman integrals to the case of massive quark lines with local operator insertions 
has been used to calculate ladder- and Benz-topologies [23,24], also leading to new types of finite 
nested sums extending the harmonic [26], generalized harmonic [27,28], and cyclotomic sums [29] 
and the associated polylogarithms, cf. Section 4. Section 5 contains the conclusions. 

2. Gluonic OMEs O^T^Cf a ) 

The calculation of all 0{nfTpCjr^) contributions to the massive OMEs has been completed with 
the computation of A ggt Q and A gq _Q at this order in [22]. In these and other computations described 

'Present analyses use the NLO corrections in x-space [5] or Mellin space implementations [6]. 
2 The O(a^) and O(a^e) contributions, with e = D — 4 the dimensional parameter, needed in the renormalization 
were given in [7, 18, 19]. 



2 



Three-Loop Contributions to the Gluonic Massive OMEs 



Johannes Bliimlein 



below we used the codes QGRAF [30], Form and tf orm [31], and color [32]. For the check of 
individual moments of the expressions derived we also used MAT AD [33]. 

In the calculation of the 0(n /P^CV^-terms large amounts of nested sums emerge. An im- 
portant step in the calculation consists in merging these individual sums to a smaller amount of 
sums containing very voluminous summands which are then solved with the codes used within 
the package Sigma [34], like EvaluateMultiSum and SumProduction, [35], written in 
mathematica [36], cf. also [37]. Moreover, this compactification shall be performed for whole 
diagrams to avoid the intermediate emergence of a larger amount of generalized harmonic sums as 
was observed in [21]. Generalized harmonic sums do not contribute in this case. 

The constant contributions to the unrenormalized massive OMEs A gg g and A gq ^Q, UgjQ, j = 
q,g, read : 

r2 J l6(N 2 +N + 2) (I 2 
Sj = C ' T *»f[- 9iN - mN+1) (j 5 * + + 3 53 + 14 & + 

16(8Ar 3 + 13jV 2 + 277y+16) 2 
+ 21(N-l)N(N+l) 2 (3?2 + Sl +Sl} 

32 (35N 4 + 91N 3 + 178/V 2 + 18(W + 70) 

21(N-l)N(N+\y Sl 

32 (U3SN 5 + 4237N 4 + 886 IN 3 + 11668N 2 + 8236N + 2276) 
+ 243(N-l)N(N+l) 4 



(3),n,7£ _ T 2 I r 1 
W - n f 1 F\ L A {N _ l){N + 2) 



APi 2 8P 2 



21N 2 {N+\) 2 1 129N 3 {N+l) 3 



+^-(N-l)(N + 2)i; 2 S l - t^V- l)(AT+2)C 3 5i + ; 



27 v /v 27' ' y ^ 729N 4 (Af+l) 4 

2P 4 56(3N 4 + 6N 3 + 13N 2 + 10N+16) 4P 5 



21N 2 {N+l) 2 ^ 21N{N+l) ' 27Ar 2 (Ar+l) 



(#-!)(# + 2) 



H2(Ar2 + jv + 2 ) 2 16P 6 
27AT 2 (jV+l) 2 1 27W 3 (Af+l) 3 1 



32P 7 i 6 (jy2 +jV + 2 )2 i 6 (jy2 + jV + 2 )2 

r^i H T.77/A7 i i\9 — W^i H TA77/.7 i i\9 — ^2^1 



8UV 4 (W+1) 4 1 3iV 2 (^+l) 2 ^ z 1 3^V 2 (A^+1) 2 



32P 8 16P 9 448(Af 2 +N + 2) 2 16P, 



243N 5 {N+\) 5 9N 3 {N+\) 3 ^ 9N 2 {N+\) 2 w 9Ar 3 (Af+l) 3 

(2.2) 



160(iV 2 +iV + 2) 2 



27AT 2 (Af+l) 2 

The corresponding 1/e terms contain contributions of the 3-loop anomalous dimensions [38] which 
we have verified by an independent calculation. Moreover a prediction made in [39] has been 
confirmed. 

The gluonic OMEs A ggt Q and A gq ^Q are needed for correct flavor matching in case of the tran- 
sition of a single heavy quark becoming light, cf. [3, 10]. Here the correct choice of the matching 
scale is of importance [40]. 
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In extending the above calculation the computation of some topological classes of diagrams 
has been automated mapping the graph to expressions involving hypergeometric p+iFp-functions. 
The £-expansion leads to nested sums being calculated using the package Sigma [34]. The au- 
tomation will cover other classes soon, requiring more involved ways to match the initial functions. 

3. Gluonic OMEs 0(T 

fCf,a) an d OMEs with massive fermion lines of two different 

masses 



All basic scalar topologies contributing to the 0(TpCp,A) of A ggi Q and A gq Q have been calculated. 



An example is given in (3.1) for the graph containing two massive triangles with m\ = ni2, 




1 (27Af 2 -5Af+16) f2N 
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Here we chose a minimal representation of sums being pairwise transcendent, which is proven by 
Sigma [34]. In these and similar diagrams sums of the kind 



N Ai 



1=1 
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OW3 \ 



(3.2) 



5^72A 5l(0 = Jo dx ^^i U 0m( x )-KA x )-KJx)-2H2)KM (3-3) 



emerge, which can be represented as Mellin transforms of iterated integrals, extending the usual 
harmonic polylogarithms (HPLs) [41], where h;(jc) = 1,H* ~(x) = dyf b (y)Hi(y),H;(x) = ln(l- 
x), and / W3 (x) = 1/ (x^l —x), cf. [42]. The relative transcendence of the respective HPLs is proven 
using differential field methods [43], cf. also [44]. In exceptional cases they may be obtained in 



terms of HPLs with root arguments. Eq. (3.1) obviously is recurrent in N. The asymptotic expan- 



sions of (3.2, 3.3) are given by 
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with N = Ne rE and Je the Euler-Mascheroni number. The poles in the complex plane are situated at 



non-positive integers and half-integers and ( p.l| ) is thus a meromorphic function. The half-integer 
poles emerge algebraically and from structures like 

'2N\ 1 _ r(iV+l/2) 



(3.6) 



N J 4 N y/nT(N+l) ' 
With these properties it can be defined in the complex A^-plane. By similar arguments the analytic 
continuation for the whole TpCpA -contribution to A gg _Q is obtained. 

First results on massive OMEs with two fermion lines for m\ / 1112 have been reported in 
[12, 13] for the moments ,/V = 2,4,6 for Aq s . The calculation has been performed by mapping the 
OMEs to tadpoles which were computed using qexp [45]. The renormalization of these matrix 
elements generalizes the case of n/ massless and one massive fermion [3] and is given in Ref. [14]. 
Since ml /ml ~ 1/10, charm cannot be treated as massless at the scale m|, which makes the use 
of the variable flavor scheme, see [3], built on this assumption, veiy problematic in case of these 
contributions. On the other hand, the fixed flavor number scheme which can be used in precision 
deep-inelastic world data analyses, cf. [1], can naturally accommodate these terms. The corre- 
sponding calculation in case of general values of N is underway [46]. 

4. Ladder and Benz Topologies 

In the following we discuss higher topologies which contribute to the massive OMEs A ; y at 3-loop 
order. These are ladder-, Benz-, and crossed box topologies with the respective local operator 
insertions (of up to four lines). The basic ladder topologies have been calculated in Ref. [23] up to 
3-leg operator insertions. Here we consider the case of only one massive fermion line. As has been 
described in [47] the Feynman diagrams can be represented as multiply nested sums. The ladder 
diagrams with six massive propagators have representations in terms of the Appell function F\ [48]. 
Most of the integrals can be solved using Sigma [34], including the pole structure in e. This is 
presently more involved in case of operator insertions with more than two legs at six massive lines. 

For the non-divergent graphs the extension of the method [25] to local operators and massive 
lines allows the calculation. We consider the diagram shown in Fig. 1 . 
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Figure 1: 3-loop ladder diagram containing a 3-vertex local operator insertion. 
The local operator insertion can be resummed introducing a subsidiary parameter x 

1 4a 1 4b ^ £-1 4a 1 \b ~ L T —T 

7=0 N=0j=0 N=0 i4 « 14b 



1 



T 4a — 



1 1 



l-xT 4a l-xT 4b 



(1 -xT 4a )(l -xT 4h ) 



(4.1) 
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The integral may then be performed using the method [25] and is expressed in terms of hyperloga- 
rithms Lg(x). They obey the relations 

r x dv 1 
UA X )= -tL 3 (j), Le(y) = l; L ,..., x) = -ln n (x , a u b e R . (4.2) 
Jo y — b n\ 

n 

As an intermediary result one obtains 



l+x lx-l 3(1 -x) l-lx + x 2 l-2x 2 
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xr 

(1 +x) 1 ^ 
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2x 3 



[3Li/2,o,o,i — Li/2,0,1,1 ~~ 3L 1 / 2 ,i,o,i +Li/2,i,i,i] 
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5 5 T 3(3 +x) 

M),o,o,i,i — XT — r- — 7^0,0,1,0,1 + — — — v-m),o,i,i,i 



-1+jc 2(-l+x) u ' u,1 ' u ' 1 2(-l+x)x 

(-l + 2x 2 ) 

2^ [3L 01 /2,o,o,i +Lo,i/2,o,i,i +3L j/ 2 ,i,o,i — L ,i/2,i,i,iJ 

3(l-3x 2 + 3x 3 ) 8- 14x + 5x 2 + 3x 3 

+ 2(-l+x)x 3 L °4A0,i+ 2( _ 1+x)x 3 Lo.i.o.1,1 

8-15x + 3x 2 3(-3 + 2x) -6 + 3x + 5x 2 
+ 2(-l+x)x 3 L o.i,i,o,i 2^—1^,1,1,1,1 + ~, L 1A o,i,, 

2(-l+x) 4-2x + 5x 2 -4 + 6x + 3x 2 
H ^3 Li, 1,1,0,1 H 2^ Li,o,i,o,i 2^3 Li,o,i,i,i 

(-1+2X 2 ) 3(-l+x)(4 + 3x) T 

H 2^ L !> 1/2,0,0,1 — M, 1/2,0,1,1 J 2^3 M, 1,0,0,1 

R+ 2 * 2 ) r T , T 1 (~l+x)(5 + 3x) 

2^3 L L 1, 1/2,1,0,1 + M,l/2,l,l,lJ 2^3 L l, 1,0,1,1 • ( 4 -->) 

Here iterated integrals over the alphabet {0,1,-1,1/2} contribute. Now the Mh Taylor coefficient 
has to be obtained for U{x), which is possible using the package Harmonic Sums [28] : 

^ = 2(1+ AT) 5 (2 + AO 5 (3+ AO 5 + (1+/V) 2 (2 + /V) 2 (3+/V) 2 ^ 3 



l) N {<S5 + 1^ + 56^ + 1^ +N A ) {-2A-5N + 2N 1 ) 3 

2(1+/V) 2 (2 + /V) 2 (3+/V) 2 + 12(2 + /V) 2 (3+/V) 2>Sl 



1 _o 1 



2 2 

2(1 +N)(2 + N){3+N) Sl + {2 + N)(3+N) SlS2 
3 14 + 63 IN + 578N 2 + 288N 3 + 6SN 4 + 5N 5 2 3 
+ 4(1+/V) 3 (2 + /V) 2 (3+/V) 2 1_ 2 5 
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(399 + 2069/V + 211 AN 2 + 151(W 3 + 349/V 4 + 27/V 5 ) 

S3 — 2S-2.-3 



6(1+/V) 2 (2+A0 2 (3+N) 2 

{-\) N (65 + WIN + 56N 2 + 13/V 3 + /V 4 ) 
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4(1+AT) 3 (2 + A0 2 (3+A0 2 2 ~ 2 

3 (99 + 225/V+19(W 2 + 65/V 3 + 7/V 4 ) 
--S3S2 - 2S 2A S 2 + 2(1 +^2(2+^)2(3+^) ^ 

1 (11 + 5AQ 

(1+/V) 4 (2 + /V) 4 (3+A0 4 1 (1+N)(2 + N)(3+Ny 3 1 
(470 + 1015N + 996/V 2 + 447/V 3 + 96N 4 + 8/V 5 ) 

4(1+/V) 2 (2 + A0 2 (3+/V) 2 5 ' 51 " 52,3 

(53 + 29/V) 3(3 + 2jV) 

+ 2(l+iV)(2 + iV)(3 + iV) 3 1 (1+N)(2 + N)(3+N) 1 U 
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2 /v (-28-25Af-4Ar 2 + Ar 3 ) /l 
+55 2A1 +65 3 ,i,,+ (1 + A0 2( 2 + iV )(3+ A0 2 5 M,i l^' 1 ' 1 

( 5 +#) c ^ 1 ^ ( 5 +^) c 1 



(1 +AT)(3 +AT) ^ V 2 ' 2 ' V " 2(l+A0(3+iV) 51 - 1 ' 1 ' 1 2' *' \J ' ^ 

Here P; denote polynomials, cf. [23]. The final expression contains individual terms which grow 
°c 2 W for /V — > 00. However, this singularity cancels in t^(N). As an extension of the usual harmonic 
sums [26] also generalized sums [27,28] contribute with weights x, £ {1, 1/2,2}. 

Let us now turn to the graph shown in Fig. 2. It consists out of two contributions with respect 
to the operator insertion, cf. [3], which may be viewed as being obtained by contraction of the 
central lines of the ladder graph, resp. the graph of the crossed box, with central operator insertion. 
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Figure 2: 3-loop ladder diagram containing a 4-vertex local operator insertion. 

While the calculation of the former graph is rather straightforward, in the latter case root-letters ap- 
peal - in the integration formalism using the method of hyperlogarithms. Here it is possible, however, 
to move the corresponding root-expressions completely into the argument of the hyperlogarithms. 
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In the next step we would like to determine the Nth Taylor coefficient of the expression obtained. 
A first possibility consists in calculating a large number of Mellin-moments for the expression in 
an efficient way, which needs the use of Form [31] beyond the representations in Maple [49]. In 
the present case about 1500 moments have been calculated. The method of guessing [50] allows 
to derive a corresponding difference equation, for which we needed ~ 700 moments in the present 
case. This equation can now be solved using Sigma [34] and the Mh moment is obtained. In 
an earlier investigation we have determined all 3-loop anomalous dimensions and massless Wilson 
coefficients from their moments in this way, which required 5114 moments, [51]. Although not 
expected to fail with a significant probability, and having a large verification space with ~ 1500 
moments available at ~ 700 needed, still the Mh moment shall be derived having a proof certifi- 
cate. A more involved calculation using HarmonicSums [28] and Sigma [34] provides this. 
Here, difference equations of up to order = 16 and degree d = 1 08 have to be solved. The results 
are presented in [24]. 

The emergence of root-expressions in x implies in the present result quite a series of new sums, 
which are of the nested binomial- and inverse-binomial type and also contain generalized harmonic 
sums. We have translated the result into x-space, where new iterated integrals with various root- 
type letters were obtained, extending those having appeared in [52]. Using the methods of Ref. [43] 
the relative transcendence of the different functions can be checked and we derived a corresponding 
basis for these functions, cf. [42]. In Af-space the basis representation, including the new sums, is 
derived using Sigma, [34]. 

The method of hyperlogarithms also allows the calculation of non-singular Benz-graphs, which 
have representations in terms of harmonic sums and generalized harmonic sums, cf. [24, 53]. 

5. Conclusions 

We reported on recent progress in calculating the asymptotic heavy flavor Wilson coefficients con- 
tributing to the deep-inelastic structure function F2(x,Q 2 ) at 3-loop order. As a first class all con- 
tributions of 0(tifTpCA.F) have been completed and basis integrals for the class 0(TpCa,f) were 
obtained. Furthermore, the automatic calculation of related topology classes for other color fac- 
tors started. We extended the former analysis to the case in which two heavy quarks of different 
mass contribute and obtained a series of Mellin moments. These contributions are no longer in 
accordance with the variable flavor number scheme, since charm does not decouple at the mass 
scale of the bottom quark. In case of two massive fermion lines with equal mass new sums and 
iterated integrals appear beyond the usual harmonic sums and polylogarithms, also leading to more 
singularities in the complex plane. We have also calculated the scalar graphs contributing to ladder- 
topologies for up to six massive fermion lines, including graphs with local 4-leg operators. Here 
the results are given in terms of special classes of generalized sums, which individually may even 
diverge exponentially for N — > °°. This divergence is canceled between different terms contributing. 
In case of the graph shown in Fig. 2 a larger amount of nested binomial and inverse binomial sums 
weighted with harmonic sums and their generalization, also including cyclotomic sums, emerge, 
extending the known alphabets both for the sums and the associated iterated integrals. This all is 
invisible at the level of Mellin moments since the corresponding expressions are given in terms of 
rationals and single ^-values. It appears that in the presence of a single mass already at 3-loop 
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order rich new structures are contributing in the single differential case being characterized either 
by the Mellin variable N or the momentum fraction x. 

Acknowledgment. We would like to thank F. Brown for discussions. The Feynman diagrams 
have been drawn using Axodraw [54]. This work has been supported in part by DFG Sonder- 
forschungsbereich Transregio 9, Computergestiitzte Theoretische Teilchenphysik, by the Austrian 
Science Fund (FWF) grant P203477-N18, by the EU Network LHCPHENOnet PITN-G A-20 1 0- 
264564, and ERC Starting Grant PAGAP FP7-257638. We thank WolframResearch for their kind 
cooperation within LHCPHENOnet. 

References 

[1] S. Alekhin, J. Bliimlein and S. Moch, Parton Distribution Functions and Benchmark Cross Sections 
atNNLO, Phys. Rev. D 86 (2012) 054009 [arXiv: 1202.2281 [hep-ph]]. 

[2] S. Bethke et ai, Workshop on Precision Measurements ofa s , arXiv:l 1 10.0016 [hep-ph]. 

[3] I. Bierenbaum, J. Bliimlein and S. Klein, Mellin Moments of the O(a^) Heavy Flavor Contributions 
to unpolarized Deep-Inelastic Scattering at Q 2 3> m 2 and Anomalous Dimensions, Nucl. Phys. B 820 
(2009)417, [arXiv:0904.3563 [hep-ph]]. 

[4] J. Bliimlein, S. Klein and B. Todtli, 0(a 2 ) and 0(0$) Heavy Flavor Contributions to Transversity at 
Q 2 > m 2 , Phys. Rev. D 80 (2009) 094010 [arXiv:0909.1547 [hep-ph]]. 

[5] E. Laenen, S. Riemersma, J. Smith and W. L. van Neerven, Complete 0(a s ) corrections to heavy 
flavor structure functions in electroproduction Nucl. Phys. B 392 (1993) 162; 229; 
S. Riemersma, J. Smith and W. L. van Neerven, Rates for inclusive deep inelastic electroproduction of 
charm quarks at HERA, Phys. Lett. B 347 (1995) 143 [hep-ph/941 1431]. 

[6] S. I. Alekhin and J. Bliimlein, Mellin representation for the heavy flavor contributions to deep 
inelastic structure functions, Phys. Lett. B 594 (2004) 299 [hep-ph/0404034]. 

[7] M. Buza, Y. Matiounine, J. Smith, R. Migneron and W. L. van Neerven, Heavy quark coefficient 
functions at asymptotic values Q 2 > m 2 , Nucl. Phys. B 472 (1996) 611 [hep-ph/9601302]. 

[8] J. A. M. Vermaseren, A. Vogt and S. Moch, The Third-order QCD corrections to deep-inelastic 
scattering by photon exchange, Nucl. Phys. B 724 (2005) 3 [hep-ph/0504242]. 

[9] S. Alekhin, J. Bliimlein, K. Daum, K. Lipka and S. Moch, Precise charm-quark mass from 
deep-inelastic scattering, arXiv: 1212.2355 [hep-ph]. 

[10] M. Buza, Y. Matiounine, J. Smith and W. L. van Neerven, Charm electroproduction viewed in the 
variable flavor number scheme versus fixed order perturbation theory, Eur. Phys. J. C 1 (1998) 301 
[hep-ph/96 12398]. 

[11] I. Bierenbaum, J. Bliimlein and S. Klein, The Gluonic Operator Matrix Elements at 0(a 2 ) for DIS 
Heavy Flavor Production, Phys. Lett. B 672 (2009) 401 [arXiv:0901.0669 [hep-ph]]. 

[12] J. Ablinger, J. Bliimlein, S. Klein, C. Schneider and F. WiBbrock, 3-Loop Heavy Flavor Corrections 
to DIS with two Massive Fermion Lines, arXiv: 1 106.5937 [hep-ph]. 

[13] J. Ablinger, J. Bliimlein, A. Hasselhuhn, S. Klein, C. Schneider and F. WiBbrock, New Heavy Flavor 
Contributions to the DIS Structure Function F 2 (x, Q 2 ) at 0(af), PoS (RADCOR201 1) 031 
arXiv: 1202.2700 [hep-ph]. 



9 



Three-Loop Contributions to the Gluonic Massive OMEs 



Johannes Bliimlein 



[14] J. Bliimlein and F. WiBbrock, in preparation. 

[15] I. Bierenbaum, J. Bliimlein and S. Klein, Logarithmic 0(af) contributions to the DIS Heavy Flavor 
Wilson Coefficients at Q 2 > m 2 , PoS DIS 2010 (2010) 148 [arXiv: 1008.0792 [hep-ph]]. 

[16] I. Bierenbaum, J. Bliimlein, S. Klein, and F. WiBbrock, in preparation. 

[17] H. Kawamura, N. A. Lo Presti, S. Moch and A. Vogt, On the next-to-next-to-leading order QCD 
corrections to heavy-quark production in deep-inelastic scattering, Nucl. Phys. B 864 (2012) 399 
[arXiv: 1205.5727 [hep-ph]]. 

[18] I. Bierenbaum, J. Bliimlein and S. Klein, Two-Loop Massive Operator Matrix Elements and 

Unpolarized Heavy Flavor Production at Asymptotic Values Q 2 ^> m 2 , Nucl. Phys. B 780 (2007) 40 
[hep-ph/0703285]. 

[19] I. Bierenbaum, J. Bliimlein, S. Klein and C. Schneider, Two-Loop Massive Operator Matrix Elements 
for Unpolarized Heavy Flavor Production to 0{e), Nucl. Phys. B 803 (2008) 1 [arXiv:0803.0273 
[hep-ph]]. 

[20] J. Bliimlein, A. De Freitas, W. L. van Neerven and S. Klein, The Longitudinal Heavy Quark Structure 
Function F® Q in the Region Q 2 > m 2 at O(o? s ), Nucl. Phys. B 755 (2006) 272 [hep-ph/0608024]. 

[21] J. Ablinger, J. Bliimlein, S. Klein, C. Schneider and F. WiBbrock, The O(a^) Massive Operator 

Matrix Elements of 0{N f) for the Structure Function F2(x,Q 2 ) and Transversity, Nucl. Phys. B 844 
(2011) 26 [arXiv: 1008.3347 [hep-ph]]. 

[22] J. Bliimlein, A. Hasselhuhn, S. Klein and C. Schneider, The O(a^nfTpCA.F) Contributions to the 
Gluonic Massive Operator Matrix Elements, Nucl. Phys. B 866 (2013) 196 [arXiv: 1205.4184 
[hep-ph]]. 

[23] J. Ablinger, J. Bliimlein, A. Hasselhuhn, S. Klein, C. Schneider and F. WiBbrock, Massive 3-loop 
Ladder Diagrams for Quarkonic Local Operator Matrix Elements, Nucl. Phys. B 864 (2012) 52 
[arXiv: 1206.2252 [hep-ph]]. 

[24] J. Ablinger, J. Bliimlein, C. Schneider, and F. WiBbrock, in prepration. 

[25] F. Brown, The Massless higher-loop two-point function, Commun. Math. Phys. 287 (2009) 925 
[arXiv:0804.1660 [math.AG]]. 

[26] J. A. M. Vermaseren, Harmonic sums, Mellin transforms and integrals, Int. J. Mod. Phys. A 14 (1999) 
2037 [hep-ph/9806280]; 

J. Bliimlein and S. Kurth, Harmonic sums and Mellin transforms up to two loop order, Phys. Rev. D 
60 (1999) 014018 [hep-ph/9810241]. 

[27] S. Moch, P. Uwer and S. Weinzierl, Nested sums, expansion of transcendental functions and 
multiscale multiloop integrals, J. Math. Phys. 43 (2002) 3363 [hep-ph/01 10083]; 

[28] J. Ablinger, J. Bliimlein and C. Schneider, Analytic and Algorithmic Aspects of Generalized 
Harmonic Sums and Polylogarithms, DESY 12-210; 

J. Ablinger, Computer Algebra Algorithms for Special Functions in Particle Physics, J. Kepler 
University Linz. PhD Thesis. April 2012; 

http : / /www . rise . jku . at /re search /combinat / software/ 

[29] J. Ablinger, J. Bliimlein and C. Schneider, Harmonic Sums and Polylogarithms Generated by 
Cyclotomic Polynomials, J. Math. Phys. 52 (2011) 102301 [arXiv: 1105.6063 [math-ph]]. 

[30] P. Nogueira, Automatic Feynman graph generation, J. Comput. Phys. 105 (1993) 279. 



10 



Three-Loop Contributions to the Gluonic Massive OMEs 



Johannes Bliimlein 



[31] J. A. M. Vermaseren, New features of FORM, math-ph/00 10025; 

M. Tentyukov and J. A. M. Vermaseren, The Multithreaded version of FORM, Comput. Phys. 
Commun. 181 (2010) 1419 [hep-ph/0702279 [HEP-PH]]. 

[32] T. van Ritbergen, A. N. Schellekens and J. A. M. Vermaseren, Group theory factors for Feynman 
diagrams, Int. J. Mod. Phys. A 14 (1999) 41 [arXiv:hep-ph/9802376]. 

[33] M. Steinhauser, MATAD: A Program package for the computation ofMAssive TADpoles, Comput. 
Phys. Commun. 134 (2001) 335 [hep-ph/0009029]. 

[34] C. Schneider, A refined difference field theory for symbolic summation J. Symbolic Comput. 43 (2008) 
611 [arXiv:0808. 2543]; Product representations in YlL-fields Ann. Comb. 9 (2005) 75; Solving 
parameterized linear difference equations in terms of indefinite nested sums and products, J. Differ. 
Equations Appl. 11 (2005) 799; Parameterized Telescoping Proves Algebraic Independence of Sums 
Ann. Comb. 14 (2010) 533 [arXiv:0808.2596]; A Symbolic Summation Approach to Find Optimal 
Nested Sum Representations. In A. Carey, D. Ellwood, S. Paycha, and S. Rosenberg, editors, Motives, 
Quantum Field Theory, and Pseudodifferential Operators, volume 12 of Clay Mathematics 
Proceedings 285, Amer. Math. Soc, 2010 arXiv:0808.2543; Symbolic summation assists 
combinatorics Sem. Lothar. Combin., 56 (2007) 1, Article B56b; Multisummation in Difference 
Fields, Habilitationsschrift JKU Linz (2007) and references therein; 
J. Ablinger, J. Bliimlein, S. Klein and C. Schneider, Modem Summation Methods and the 
Computation of 2- and 3-loop Feynman Diagrams, Nucl. Phys. Proc. Suppl. 205-206 (2010) 110 
[arXiv: 1006.4797 [math-ph]]. 

[35] C. Schneider, in preparation. 

[36] http : / / www . wolfram . com/ mathemat ica/ 

[37] J. Bliimlein, A. Hasselhuhn and C. Schneider, Evaluation of Multi-Sums for Large Scale Problems, 
arXiv: 1202.4303 [math-ph]; in: Proceedings of 10th International Symposium on Radiative 
Corrections, PoS(RADCOR201 1)32, (2012); 

J. Ablinger, S. Bliimlein, M. Round and C. Schneider, Advanced Computer Algebra Algorithms for 
the Expansion of Feynman Integrals, arXiv:1210.1685 [cs.SC]. 

[38] A. Vogt, S. Moch and J. A. M. Vermaseren, The Three-loop splitting functions in QCD: The Singlet 
case, Nucl. Phys. B 691 (2004) 129 [hep-ph/04041 1 1]. 

[39] J. E Bennett and J. A. Gracey, Determination of the anomalous dimension of gluonic operators in 
deep inelastic scattering at 0(1 /N f ), Nucl. Phys. B 517 (1998) 241 [hep-ph/97 10364]. 

[40] J. Bliimlein and W. L. van Neerven, Heavy flavor contributions to the deep inelastic scattering sum 
rules, Phys. Lett. B 450 (1999)417 [hep-ph/981 1351]. 

[41] E. Remiddi and J. A. M. Vermaseren, Harmonic polylogarithms, Int. J. Mod. Phys. A 15 (2000) 725 
[hep-ph/9905237]. 

[42] J. Ablinger, J. Bliimlein, C. Raab and C. Schneider, in preparation. 

[43] M. Bronstein, Symbolic Integration I - Transcendental Functions, (Springer, Heidelberg, 1997). 

[44] C. Raab, Definite Integration in Differential Fields, PhD Tesis, August 2012, Johannes Kepler 
University Linz, Austria. 

[45] R. Harlander, T. Seidensticker, M. Steinhauser, Complete corrections ofO(aa s ) to the decay of the 
Z-boson into bottom quarks, Phys. Lett. B426 (1998) 125, [hep-ph/9712228]; 
T. Seidensticker, Automatic application of successive asymptotic expansions of Feynman diagrams, 
[hep-ph/9905298]. 



11 



Three-Loop Contributions to the Gluonic Massive OMEs 



Johannes Bliimlein 



[46] J. Bliimlein, M. Round, C. Schneider, F. WiBbrock, in preparation. 

[47] J. Bliimlein, Structural Relations of Harmonic Sums and Mellin Transforms up to Weight w — 5, 
Comput. Phys. Commun. 180 (2009) 2218 [arXiv:0901.3106 [hep-ph]]; 
J. Bliimlein, S. Klein, C. Schneider and F. Stan, A Symbolic Summation Approach to Feynman 
Integral Calculus, J. Symbolic Comput. 47 (2012) 1267-1289 [arXiv: 101 1.2656 [cs.SC]]. 

[48] P. Appell, Sur Les Fonctions Hypergeometriques de Plusieurs Variables, (Gauthier- Villars, Paris, 
1925); 

P. Appell and J. Kampe de Feriet, Fonctions Hypergeometriques; Polynomes d'Hermite, 
(Gauthier- Villars, Paris, 1926); 

W.N. Bailey, Generalized Hypergeometric Series, (Cambridge University Press, Cambridge, 1935); 
A. Erdelyi (ed.), Higher Transcendental Functions, Bateman Manuscript Project, Vol. I, 
(McGraw-Hill, New York, 1953); 

L.J. Slater, Generalized Hypergeometric Functions, (Cambridge University Press, Cambridge, 1966); 
H. Exton, Multiple Hypergeometric Functions and Applications, (Ellis Horwood Limited, Chichester, 
1976); 

H. Exton, Handbook of Hypergeometric Integrals, (Ellis Horwood Limited, Chichester, 1978). 

[49] http : / / www . maplesof t . com/products /maple/ 

[50] M. Kauers, Guessing Handbook, Technical Report RISC 09-07 (2009), JKU Linz. 

[51] J. Bliimlein, M. Kauers, S. Klein and C. Schneider, Determining the closed forms of the O(a^) 

anomalous dimensions and Wilson coefficients from Mellin moments by means of computer algebra, 
Comput. Phys. Commun. 180 (2009) 2143 [arXiv:0902.4091 [hep-ph]]. 

[52] U. Aglietti and R. Bonciani, Master integrals with 2 and 3 massive propagators for the 2 loop 
electroweak form-factor - planar case, Nucl. Phys. B 698 (2004) 277 [hep-ph/0401 193]. 

[53] J. Ablinger, J. Bliimlein, A. De Freitas, A. Hasselhuhn, S. Klein, C. Schneider and F. WiBbrock, New 
Results on the 3-Loop Heavy Flavor Wilson Coefficients in Deep-Inelastic Scattering, 
arXiv: 1212.5950 [hep-ph], PoS(ICHEP2012)270. 

[54] J. A. M. Vermaseren, Axodraw, Comput. Phys. Commun. 83 (1994) 45. 



12 



